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Abstract. The El Nino Southern Oscillation is the domi- 
nant year-to-year mode of global climate variability. El Nino 
effects on terrestrial carbon cycling are mediated by associ- 
ated climate anomalies, primarily drought, influencing fire 
emissions and biotic net ecosystem exchange (NEE). Here 
we evaluate whether El Nino produces a consistent response 
from the global carbon cycle. We apply a novel bottom- 
up approach to estimating global NEE anomalies based on 
FLUXNET data using land cover maps and weather reanal- 
ysis. We analyze 13 years (1997-2009) of globally gridded 
observational NEE anomalies derived from eddy covariance 
flux data, remotely-sensed fire emissions at the monthly time 
step, and NEE estimated from an atmospheric transport in- 
version. We evaluate the overall consistency of biospheric 
response to El Nino and, more generally, the link between 
global CO 2 flux anomalies and El Nino-induced drought. 
Our findings, which are robust relative to uncertainty in both 
methods and time-lags in response, indicate that each event 
has a different spatial signature with only limited spatial co- 
herence in Amazonia, Australia and southern Africa. For 
most regions, the sign of response changed across El Nino 
events. Biotic NEE anomalies, across 5 El Nino events, 
ranged from -1.34 to +0.98PgCyr _1 , whereas fire emis- 
sions anomalies were generally smaller in magnitude (rang- 
ing from -0.49 to +0.53PgCyr _1 ). Overall drought does 
not appear to impose consistent terrestrial CO 2 flux anoma- 
lies during El Ninos, finding large variation in globally in- 
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tegrated responses from -1.15 to +0.49 Pg C yr“ 1 . Despite 
the significant correlation between the CO 2 flux and El Nino 
indices, we find that El Nino events have, when globally inte- 
grated, both enhanced and weakened terrestrial sink strength, 
with no consistent response across events. 


1 Introduction 

The El Nino/Southern Oscillation (ENSO), a global atmo- 

spheric circulation feature linked to atmospheric pressure 
patterns and sea surface temperature anomalies in the trop- 
ical Pacific (McPhaden et al., 2006), is the dominant mode 
of global year-to-year climate variation (Buermann et ah, 
2003; Trenberth et ah, 2007). Although ENSO originates in 
the equatorial east Pacific ocean it influences terrestrial car- 
bon cycling globally through teleconnections to water bal- 

ance (McPhaden et ah, 2006; Ropelewski and Halpert, 1986, 
1987; Reichenau and Esser, 2003; Williams and Hanan, 
2010; Woodward et ah, 2008), temperature (McPhaden et 

ah, 2006; Nagai et ah, 2007), and fire emissions (Page et ah, 
2008; van der Werf et ah, 2004, 2006). The net effect of 
ENSO on terrestrial carbon cycling is assumed to be linked 
to its phase, with the warmer El Nino phase commonly as- 
sociated with a decrease in terrestrial uptake of CO 2 (e.g., 
Gurney et ah, 2008; Jones et ah, 2001; McGuire et ah, 2001; 
Rodenbeck et ah, 2003; Qian et ah, 2008); although Bosquet 
et ah (2000) showed both anomalous CO 2 uptake and efflux 
during El Nino years. 
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The tropics, the most dominant region in terms of inter- 
annual variability of global CO 2 flux (Rayner et al., 2008), 
plays an important role in modulating the effect of El Nino 
on carbon cycling. El Nino events have been reported to 
have distinct spatial signatures with coherent responses, e.g., 
warmer and drier conditions in Tropical Asia during the dry 
season (Woodward et al., 2008), droughts and reduced net 
carbon uptake in Amazonia (Aragao et al., 2007), and en- 
hanced fire emissions, especially in insular Southeast Asia 
(van der Werf et al., 2006). However, observations and 
model-based studies do not agree on whether El Nino and 
associated drought elicits a regularized response in terrestrial 
sink strength. In Amazonia terrestrial biosphere models have 
been shown to invert the annual carbon flux cycle relative to 
observed data (Saleska et al., 2003), i.e., models incorrectly 
simulated net uptake during the wet season and net efflux 
during the dry season (cf. Aragao et al., 2007). In a sub- 
sequent biophysical model study Baker et al. (2008) identi- 
fied increased photosynthesis (in response to enhanced light 
levels) and decreased respiration (in response to surface soil 
desiccatation) during the dry season as necessary mechanis- 
tic concepts to match observed seasonality. 

At the footprint scale Schwalm et al. (2010a) observed 
a net increase in CO 2 uptake coincident with drought in 
tropical evergreen broadleaf forests. In contrast, Keller et 
al. (2004) showed that the interplay between the wet and dry 
seasons relative to CO 2 uptake/efflux varied across the Ama- 
zon Basin. On coarser spatial scales enhanced uptake coinci- 
dent with drought has been documented using remote sens- 
ing (Huete et al., 2006; Saleska et al., 2007) but this result 
is disputed due to data artifacts (Samanta et al., 2010). The 
overall ambiguity in response is confounded by sparse mon- 
itoring networks in the tropics and large interannual variabil- 
ity (Bosquet et al., 2000) relative to a small source magnitude 
with large uncertainties (Gurney et al., 2002; Stephens et al., 
2007). 

Beyond net carbon uptake by ecosystems El Nino also 
influences carbon emissions by fires (e.g., van der Werf et 
al., 2004), though the relative contribution of these two fac- 
tors remains uncertain and poorly constrained in time and 
space (Canadell et al., 2007b). In general, studies that exam- 
ine the effect of El Nino on carbon cycling rely on ecosys- 
tem process models and atmospheric inversions and empha- 
size either fire emissions (Fuller and Murphy, 2006; Langner 
and Siegert, 2009; Page et al., 2008; Patra et al., 2005b; 
Rodenbeck et al., 2003; van der Werf et al., 2004) or bi- 
otic NEE (Adams and Piovesan, 2005; Hasimoto et al., 2004; 
Reichenau et al., 2003; Woodward et al., 2008; Zeng et al., 
2005). 

While ecosystem process models explicitly incorporate 
biophysical and/or biogeochemical mechanisms they per- 
form poorly when simulating observed NEE both glob- 
ally (Schwalm et al., 2010b) and in the tropics (Saleska 
et al., 2003). Inversion approaches identify large-scale net 
sources/sinks and, despite substantial variation in absolute 


magnitude, typically report similar interannual variability 
(Gurney et al., 2008). Such consistency follows from consid- 
erable overlap in the underlying observations of atmospheric 
CO 2 concentration. The network of CO 2 observations is 
sparse (Rayner et al., 1999), making it difficult for inversions 
to resolve regional carbon fluxes, particularly in the tropics 
(Gurney et al., 2003) which dominate interannual variability 
in global CO 2 flux (Rayner et al., 2008). 

FLUXNET offers an alternative to ecosystem process 
models or atmospheric inversions to evaluate the influence 
of El Nino events on the global carbon cycle. FLUXNET 
is a network of globally distributed eddy covariance towers 
with dense data streams of observed in situ CO 2 exchange 
and ancillary variables (Baldocchi, 2008). The extended spa- 
tiotemporal coverage of FLUXNET allows for the climate 
sensitivity of terrestrial CO 2 flux to be scaled up to estimate 
global flux using various data-driven approaches (Beer et al., 
2010; Jung et al., 2009; Schwalm et al., 2010a; Xiao et al., 
2008). 

Here we investigate whether global CO 2 flux anomalies 
have a clear spatial signature and consistent response during 
El Nino events. We calculate NEE response to El Nino events 
from 1997-2009 using a novel upscaling approach based 
on micrometeorological observations of CO 2 exchange, plus 
fire emissions from a forward biogeochemical model and 
remotely-sensed estimates of area burned, fire activity, and 
plant productivity. Total biospheric response to El Nino 
events is based on monthly anomalies: NEE sensitivities to 
drought derived from FLUXNET (Schwalm et al., 2010a) are 
spatially and temporally scaled and then combined with fire 
emissions anomalies. We generate globally gridded monthly 
anomaly maps of NEE and fire emissions components of ter- 
restrial carbon sink response. Our objective is to evaluate 
spatial and temporal coherence in terrestrial carbon flux re- 
sponse during El Nino events and to quantify the degree to 
which hydroclimatic anomalies during El Nino events influ- 
ence the terrestrial carbon cycle. 

2 Methods 

Our method isolates drought-related El Nino effects on biotic 
NEE. We also account for fire as a carbon source as fire emis- 
sions have been shown to increase during El Nino events. We 
estimated these 2 components using independent sources of 
data and added them to quantify total biospheric response to 
El Nino on a monthly time step. We based the biotic NEE 
on scaled ecosystem sensitivities (Schwalm et al., 2010a) de- 
rived from FLUXNET data (Baldocchi, 2008). These sen- 
sitivities quantified the change in NEE for a unit change in 
relative water deficit anomaly where water deficit, calculated 
as the ratio of latent heat to available energy, represented 
drought (Schwalm et al., 2010a). NEE was expressed using 
the atmospheric convention, a positive value indicates out- 
gassing of CO 2 to the atmosphere. 


Biogeosciences, 8, 2493-2506, 2011 


www.biogeosciences. net/8/2493/2011/ 



C. R. Schwalm et al.: Global CO 2 flux anomalies 


2495 



Fig. 1. Inversion and FLUXNET measurement sites. Maps display location of FLUXNET (n = 238 sites, green circles) and inversion (n = 5 1 
sites, red squares) sites used in this study both globally (Panel A) and, at a higher spatial resolution, for Europe (Panel B). Note there are 2 
additional inversion sites with variable location that are not shown. 


To calculate monthly ecosystem sensitivities, monthly in- 
tegrals of tower-based NEE and water deficit were first nor- 
malized by month and site. That is, for each site’s monthly 
time series, the mean monthly cycle of both NEE and water 
deficit were normalized by subtracting the mean and divid- 
ing by the standard deviation, which gave, in units of stan- 
dard deviation, the departure from the mean value or rela- 
tive anomaly (z-scores). Next, both sets of monthly relative 
anomalies were grouped by climatic season and site-specific 
land cover class according to the International Geosphere- 
Biosphere Program land cover classification (IGBP, Love- 
land et al., 2001). We then regressed the relative flux 
anomaly on relative water deficit anomaly within each IGBP 
class and climatic season. Regression in normalized space 
corrected for spatial gradients in the FLUXNET data compi- 
lation (Lauenroth and Sala, 1992) and yielded a linear rela- 
tionship with the slope parameter representing flux sensitiv- 
ity, i.e., change in relative NEE anomaly for a unit change 
in relative water deficit anomaly (Schwalm et al., 2010a). 
These dimensionless sensitivities were transformed into NEE 
sensitivities (g Cm -2 month - 1 er -1 ), where <7 (standard de- 
viation) referred to a 1 unit change in relative water deficit 
anomaly (negative values indicated increased dryness), using 
the standard deviation of NEE in absolute units. The NEE 
sensitivities were derived using 5173 months of FLUXNET 
data (Baldocchi, 2008) across 938 site-years collected be- 
tween 1991 and 2006 from 238 measurement towers (11 dis- 
tinct land cover classes) distributed globally (Fig. 1) with the 
highest tower densities in North America and Europe (see 
Schwalm et al. (2010a) for details). 

Monthly NEE sensitivities to water deficit were then spa- 
tially scaled using a land cover scheme. In parallel to the 
IGBP designation for each flux tower, the IGBP land cover 
dataset includes fractional values for all possible 18 classes 
for the full global land surface (Loveland et al., 2001). The 
NEE sensitivity of each terrestrial pixel (1° x 1°) was either 
estimated from FLUXNET (vegetated types, Schwalm et al., 
2010a) or effectively 0 (non-vegetated types, e.g., urban and 
built up). We combined observed NEE sensitivities with 
fractional land cover for each pixel using a weighted sum 


across all IGBP class-specific sensitivities where fractional 
coverages served as weights. As NEE sensitivities were de- 
rived monthly by climatic season, e.g., December, January, 
February as winter in the northern and summer in the south- 
ern hemisphere, this resulted in 4 distinct sensitivity maps 
(Schwalm et al., 2010a), whereby the assigned sensitivities 
for each month within a climatic season were identical. 

These pixel-specific sensitivities were subsequently tem- 
porally scaled using monthly values of water deficit from 
the Modern Era Retrospective-analysis for Research and 
Applications (MERRA, Bosilovich et al., 2008). As with 
FLUXNET data, water deficit was calculated as the ratio of 
latent heat to available energy. Each pixel was normalized 
by month; each month across the 1997-2009 period had a 
mean value of 0 and a = 1. While MERRA data extends 
from 1979 - present the 1997-2009 period was chosen to 
match the temporal coverage of fire emissions data (see be- 
low). The normalization resulted in a monthly grid of water 
deficit relative anomalies for all 156 months in the 13-year 
record. These grids were combined with the climatic season- 
specific grids of NEE sensitivity to generate spatial fields of 
flux anomaly for each month: 

r/NEE 

Abiotic — : Z, (1) 

az 

where, for any given month and pixel. Abiotic 
(g Cm -2 month -1 ) is the NEE flux anomaly, <7 NEE / clz is 
the pixel- and climatic season-specific NEE sensitivity to 
water deficit (g C rrf 2 month - 1 a" 1 ), and z is the monthly 
pixel-specific relative anomaly of water deficit in a units. 

The principal strength of this approach is its reliance on 
observations of in situ CO 2 exchange as it responds to a sin- 
gle driver, water deficit. The derived sensitivities isolate hy- 
droclimatic effects while excluding potentially confounding 
covariates. This allows a conditional analysis of the carbon 
consequences of El Nino predicated on a robust relationship 
with drought. Alternatively, El Nino is analyzed as a drought 
phenomenon based on its large influence on water balance in 
general and precipitation in particular (e.g., Ropelewski and 
Halpert, 1986, 1987; Williams and Hanan, 2010). Previous 
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diagnostic attempts to upscale point-based eddy covariance 
data, i.e., to move from site-specific anecdotes to spatially 
explicit large-scale characterizations, have typically featured 
data mining algorithms (Beer et al., 2010; Jung et al., 2008, 
2009; Papale and Valentini, 2003; Vetter et al., 2008; Xiao et 
al., 2008; Yang et al., 2006, 2007). While such approaches 
may show high skill in matching observed or reference fields, 
as black-box algorithms they provide less insight into the un- 
derlying mechanisms, drought in this case. Our technique 
enables a direct examination of the El Nino effects on terres- 
trial carbon cycling. 

The second component of total biospheric response, fire 
emissions anomalies (Afi re , g Cm -2 month -1 ), was derived 
from monthly data in the Global Fire Emissions Database 
(GFED) version 3.1 (Randerson et al., 2007; van der Werf 
et al., 2006, 2010). For each grid cell the mean monthly 
seasonal cycle from 1997-2009 was removed with the re- 
mainder giving Afl re . Both fire and biotic anomalies shared 
the same spatial resolution (1° x 1°) and were added to esti- 
mate the total drought-related terrestrial carbon sink response 
( A total) by pixel to an El Nino event. All anomalies (Abiotic, 
Afl re , and A tota i) represent the monthly deviation from the 
mean monthly cycle across the full analysis period from 
1997-2009 where positive anomalies indicate increased CO 2 
outgassing to the atmosphere relative to the 13-year mean. 

El Nino events were determined using the Multivariate 
ENSO Index (MEI, Wolter and Timlin, 1998). The bi- 
monthly MEI values were indexed by their last month, 
e.g., December- January was indexed as January. Any given 
month was classified as an El Nino condition if its bimonthly 
MEI value lay in the upper tercile of all bimonthly MEI val- 
ues from December 1949-January 1950 through November 
2009-December 2009. El Nino events of 3 months or less 
were discarded. Internal gaps (MEI values outside the upper 
tercile) of 1 month were still counted as the same event. For 
each El Nino event, the terrestrial carbon sink anomaly for 
each pixel was summed across all its months. Global event 
totals were estimated by aggregating over the full terrestrial 
land surface. 

Fire emissions and biotic NEE response to El Nino may, 
with considerable across-study variability, exhibit a lagged 
response (e.g.. Fuller and Murphy, 2006; Qian et al., 2008, 
Wooster et al., 2011). We analyzed the effect of a poten- 
tially delayed response of carbon cycling to El Nino by es- 
timating lags (/mei) between MEI and flux anomalies using 
the highest (in absolute value) statistically significant cross- 
correlation and by calculating lagged correlations as well as 
global event totals using time-lags (flux response lags MEI) 
of up to 6 months. 

We estimated uncertainty for Abiotic, Afl re , and A to tai as 
one standard deviation about the mean value (ler). For 
Abiotic, 2 sources of uncertainty were considered: uncer- 
tainty in NEE sensitivity and uncertainty in MERRA reanal- 
ysis field. The latter was based on the dispersion of multi- 
ple overlapping retrospective analyses. We used the Multi- 


model Analysis for the Coordinated Enhanced Observing 
Period dataset (CEOP, Bosilovich et al., 2009), an ensem- 
ble of 10 global reanalysis products, including MERRA, 
from October 2002-December 2004 (ftp://agdisc.gsfc.nasa. 
gov/private/ceop/), and calculated ler across all ensemble 
members for each terrestrial pixel. These values were then 
averaged by month to create a mean annual cycle of monthly 
uncertainty used for the full hindcast period. 

For the relation between NEE and water deficit, i.e., NEE 
sensitivity, we combined the uncertainties from the slope pa- 
rameter (Cook and Weisberg, 1999) and pooled standard de- 
viation of NEE (Schwalm et al., 2010). This method pro- 
vided the mean (the value used in the reported NEE anoma- 
lies) and la for both sources which were then combined 
using Monte Carlo methods as follows: (i) For all pixels 
and all months randomly draw from a Gaussian distribu- 
tion a realization of NEE sensitivity (dNEE/dz; Eq. 1) 
and water deficit, (ii) Normalize water deficit to gener- 
ate z (Eq. 1). (iii) Calculate Abiotic for all pixels and all 
months, (iv) Repeat this process 1000 times. This gener- 
ated 1000 realizations of the 156 month analysis period with 
uncertainty of Abiotic (= Abiotic) given by ler across the full 
set. As Abiotic is based on upscaled NEE sensitivities from 
FLUXNET, we also examined network representativeness by 
cross-tabulating the number of site-years used to estimate the 
NEE ~ water deficit relationship against Koppen-Geiger cli- 
mate (Peel et al., 2007) and land cover class (Loveland et al., 
2001). We then mapped this globally to obtain a map of the 
network’s degree of representation. 

For Afl re only yearly uncertainties are available. We used 
a 20% relative error for globally integrated annual integrals: 
gflre = 0.20 Afi re (van der Werf et al., 2010). Uncertainty for 
the total flux anomaly, for annual values only, combined both 

components in quadrature: e to tai = Js fi re + e biotic- Further- 
more, we created a map of ebiotic and its component sources 
of error using a weighted average across all months with 
weights based on monthly /apar (fraction of photosyntheti- 
cally active radiation absorbed by vegetation) normals, an in- 
dependent measure of seasonal variability in vegetation pro- 
ductivity. Monthly normals of /apar were derived from the 
Global Inventory Monitoring and Modeling Study (GIMMS, 
Tucker et al., 2005). 

Lastly, we compared our bottom-up flux anomalies with 
top-down estimates from the 1997-2008 Jena CO 2 in- 
version (3.75° x 5°; run s96 v3.2; update of Rodenbeck 
et al., 2003, 2005; http://www.bgc-jena.mpg.de/~christian. 
roedenbeck/download-C02/). These were based on CO 2 
concentrations from 53 measurement sites (Fig. 1) translated 
into flux anomalies with a global atmospheric tracer transport 
model. From the inversion estimates, only the globally inte- 
grated terrestrial flux was considered (fossil fuel emissions 
and ocean exchange were excluded). Anomalies were esti- 
mated similar to GFED data: the mean monthly seasonal cy- 
cle from 1997-2008 was removed with the remainder giving 
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Fig. 2. Temporal profiles of CO 2 atmospheric growth rate and MEI 
from 1997-2009. Monthly time series show MEI (red) and CO 2 
growth rate (green). Year label is centered on June-to-August. Grey 
background indicates El Nino event. 


the total NEE anomaly. Uncertainty was quantified as the 
mean absolute deviation across an ensemble of 6 runs (up- 
date of Rodenbeck et al., 2003, 2005) aggregated to annual- 
ized event totals. 


5/5 Positive 
4/5 Positive 
Equivocal 
4/5 Negative 
5/5 Negative 

4/4 Positive 
3/4 Positive 
Equivocal 
3/4 Negative 
4/4 Negative 

Fig. 3. Terrestrial carbon sink anomaly sign across multiple El Nino 
events. Panel A: Patterns of Abiotic sign (1° x 1°) f° r 5 El Nino 
events from 1997-2009. Panel B: Patterns of Jena inversion-based 
sign of NEE anomaly (3.75° x 5°) for 4 El Nino events from 1997— 
2008. Legend numbers indicate the number of events out of all 5 
(4 for Panel B ) with the relevant sign. Non- vegetated grid cells and 
those with signs split 2 to 3 (2 to 2 for Panel B) shown in white. 
A positive sign indicates increased outgassing of CO 2 to the atmo- 
sphere. 



3 Results and discussion 

Relating global CO 2 flux anomalies to El Nino requires a 
connection between biospheric flux variability and ENSO 
activity. Atmospheric CO 2 growth rate is often linked to 
ENSO activity with higher growth rates typically associated 
with the warm El Nino phase of the ENSO cycle (e.g., Gur- 
ney et al., 2008; Patra et al., 2005a). Here we related the 
Mauna Loa CO 2 monthly time series to MEI (Fig. 2). The 
Mauna Loa data was first deseasonalized and then detrended; 
this latter step removed the secular trend associated with in- 
dustrial loading and fossil fuel emissions. We found a clear 
link between ENSO and the carbon cycle at the global scale 
from 1997-2009. Variability in correlation was linked to lag 
which represented the time interval needed for changes in 
CO 2 growth rate to reach the Mauna Loa observatory (Patra 
et al., 2005a). As lag increased, correlation (p < 0.1) gener- 
ally increased: r = 0.12, 0.41, 0.55, and 0.62 for lags of 3, 6, 
9, and 12 months respectively. 

Beyond a baseline linkage between ENSO and the carbon 
cycle, the framework developed here, especially Abiotic, re- 
quires that water deficit anomalies observed with FLUXNET 
be representative of hydroclimatic conditions coincident with 
El Nino. All FLUXNET data were collected from 1991 to 
2006 with ~90% of the 938 site-years observed between 
1999 and 2006. This compares favorably to the 1997 to 
2009 analysis period. Similarly, water deficit anomalies in 
the underlying FLUXNET dataset (range: ~ ±3 a) indicated 
that FLUXNET sites observed extreme hydroclimatic condi- 
tions with several sites spatially and temporally coincident 
with the examined El Nino events. Even the tropics, rep- 
resented by the evergreen broadleaf land cover class, ex- 


perienced both abnormally dry and wet conditions (range: 
-2.85<7 to +3.28er). Furthermore, FLUXNET sites were 
sensitive to hydroclimatic anomalies with response magni- 
tudes lowest during climatic winter (e.g., nil for savannas) to 
80 g C itT 2 month - 1 for croplands during climatic summer 
(Schwalm et al., 2010a). Given the link between ENSO and 
carbon cycling as well as the sensitivity of FLUXNET sites 
to hydroclimatic anomalies, we next evaluated spatiotempo- 
ral consistency in terrestrial carbon sink response to El Nino 
during the 13-yr analysis period. 

The sign of Abiotic across all 5 El Nino events showed spa- 
tial coherence in several regions (Fig. 3a): Abiotic > 0 oc- 
curred primarily in Australasia, southern Africa, and parts of 
Amazonia. In contrast, regions where the terrestrial carbon 
sink was enhanced (Abiotic < 0) were located in the continen- 
tal United States (excluding the Upper Great Plains), tropical 
Africa, northern and southern South America as well as more 
stippled patterns throughout Eurasia. Patterns for the sign of 
Afire (not shown) were more diffuse (both signs intermixed 
without a clear emergent spatial reference) with loose clus- 
ters of positive pixels in Indonesia, Japan, and northern South 
America as well as negative pixels (less fire emissions than 
the long-term average) across much of the Northern Hemi- 
sphere. 

These coherently signed anomalies for both Abiotic and 
Afire were however small in magnitude and occurred on less 
than ~3 % of the vegetated land surface. The 5 El Nino 
events observed with GFED and scaled FLUXNET sensitiv- 
ities consisted of 55 months, excluding internal gap months. 
Using a false discovery rate of 0.05 (Ventura et al., 2004) 
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Table 1 . Flux anomaly, uncertainty, intensity, and duration of El Nino events from 1997-2009. Positive anomaly values indicate increased 
outgassing of CO 2 to the atmosphere. 


Event Flux Anomaly (Pg C yr 1 ) 


El Nino Events 

Mean MEI Value 

Duration (months) 

A biotic 

<D 

4=1 

<1 

A total 

Jena 

April 1997-June 1998 

2.16 

15 

-0.25 ±0.04 

0.53 ±0.11 

0.28 ±0.11 

2.48 ±0.15 

April 2002-April 2003 

0.83 

13 

0.39 ±0.03 

0.10 ±0.02 

0.48 ± 0.04 

1.19 ±0.03 

May 2004-May 2005 

0.59 

13 

0.13 ±0.03 

0.05 ±0.01 

0.18 ±0.04 

-0.48 ±0.07 

June 2006-February 2007 

0.83 

9 

-1.34 ±0.02 

0.19 ±0.04 

— 1.15 ± 0.05 

-1.01 ±0.02 

June 2009-December 2009 

0.94 

7 

0.98 ±0.03 

-0.49 ±0.10 

0.49 ±0.10 

— 

Mean (all events) 

1.07 

11.4 

-0.02 ±0.86 

0.08 ±0.37 

0.06 ±0.69 

- 

Mean (overlapping events) 

1.11 

12.5 

-0.27 ±0.76 

0.22 ±0.22 

-0.05 ±0.74 

0.55 ± 1.60 


Reported uncertainties represent la. Abiotic 1 uncertainties in water deficit from MERRA and FLUXNET-derived sensitivity combined using Monte Carlo methods. Afi re : 20% 
annual relative error (van der Werf et al., 2010). A tota j : uncertainties for Abiotic an d Ag re summed in quadrature. For the Jena inversion (update of Rodenbeck et al., 2003, 2005) 
uncertainties calculated as mean absolute deviation over a 6-member ensemble. Mean event rows show annualized mean and a across El Nino events. The Jena inversion record 
ends in 2008. The 2009 El Nino event continues in 2010 but GFED version 3.1 data extends only through 2009. Overlapping events are those with Abiotic an d Afire as we ll as J ena 
inversion estimates, i.e., the first four events only. 


with all El Nino months as replicates revealed 407 Ag re pix- 
els (of a possible 12 371 vegetated pixels) with non-zero 
mean values and only 4 with |Ag re | > 10 g Cm -2 month -1 . 
For Abiotic only 5 grid cells showed a non-zero mean 
(largest magnitude = -2.1 g Cm -2 month -1 ) across all El 
Nino months. 

The magnitude of non-zero anomalies coupled with their 
scarcity is strong evidence that the responses across El Nino 
events were inconsistent and not statistically robust in time 
or space. This is further supported by Jena inversion re- 
sults (Fig. 3b). These implicitly comprise fire emissions and 
biotic components and also showed minimal across-event 
consistency in carbon sink sign anomaly. Furthermore, the 
across-event inconsistency was not scale dependent. The 
consistency of the bottom-up approach was not improved 
by coarsening these values to match the Jena inversion grid 
(Fig. 4). Nor did subsequent coarsening of both top-down 
and bottom-up values (from 3.75° x 5° to 15° x 15°) reveal 
a relationship between consistency across-events and resolu- 
tion (Fig. 4). Similarly, globally aggregated El Nino anoma- 
lies from 1997-2009 (Table 1) were highly variable. This 
ambiguity of El Nino as a control on CO 2 sources/sinks was 
also present in event-based derivative totals, i.e., the time 
derivative of Abiotic during each event showed inconsistent 
sign and no relationship with sign of total Abiotic anomaly by 
event (not shown). While Ag re was generally positive (en- 
hanced outgassing of CO 2 to the atmosphere). Abiotic was 
negative for 2 of 5 events with A to tai ranging from -1.15 to 
+0.49PgCyr -1 relative to the 2000-2006 average sink of - 
2.8PgCyr -1 (Canadell et al., 2007a). Inversion-constrained 
estimates showed more variability and ranged from -1.0 to 
+2.5PgCyr -1 . 

Event-based flux anomaly totals adjusted for potential lags 
also showed no consistent response across El Nino events. 
Across a range of lags (/mei) from 0 to 6 months, globally 


aggregated El Nino flux anomalies displayed no coherent re- 
sponse (Fig. 5). Whereas Abiotic (Fig- 5a) showed the most 
variability in event-based totals, A| irc at lags of 4 months or 
more (Fig. 5b) was always positive. El Nino acted to both 
increase and decrease terrestrial carbon sink strength at all 
lags across the 5 observed events. 

Each El Nino event was largely unique in its spatial signa- 
ture. Some grid cell clusters of Abiotic consistent response 
(Fig. 6) were present across at least 3 events and corre- 
sponded to well-established zones of influence: the dipole 
pattern in sub-Saharan Africa as well as the juxtaposition of 
regions with CO 2 loss and uptake in South America. Other 
regions of biotic response have not been previously doc- 
umented: an east-west swath from Central Siberia to the 
Iberian Peninsula showed numerous pixels with | A biotic I > 
50 g C m -2 yr - 1 but of variable sign. Lastly, India and In- 
donesia, generally assumed to be drier during El Nino events, 
were also variable in sign. For Afi re (Fig. 7) patterns were 
more consistent albeit sparser. The same regions showed a 
response in all events, especially Indonesia, with changes in 
sign apparent in the northern high latitudes. 

Variability across these 5 El Nino events was large, in con- 
trast to previous studies (Gurney et al., 2008; Jones et al., 
2001; Rodenbeck et al., 2003; Qian et al., 2008), and most 
regions showed responses that switched sign across events 
(cf. Bousquet et al., 2000). This was largely confined to the 
biotic component. On a globally integrated basis Afi re shifted 
from -0. 13 to +0. 16 Pg C month - 1 during normal conditions 
to -0.09 to +0.19PgCmonth -1 across all El Nino months. 
This highlighted a general tendency of increased fire emis- 
sions during El Nino events (Table 1). The same tendency 
was not evident for Abiotic- which exhibited more variability 
outside El Nino events than during: -0.35 to +0.38 vs. -0.29 
to +0.33 Pg C month -1 respectively. 
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Fig. 4. Patterns of anomaly sign across multiple El Nino events and spatial resolutions. Left-side: Jena inversion-based sign of NEE anomaly 
at 3.75° x 5° (Panel A), 10° x 10° (Panel B) and 15° x 15° (Panel C) resolutions. Right-side: Azotic sign at 3.75° x 5° (Panel D), 10° x 10° 
(Panel E) and 15° x 15° (Panel F) resolutions. Legend numbers indicate the number of events out of all 5 (4 for Panel A-C) with the relevant 
sign. Non-vegetated grid cells and those with signs split 2 to 3 (2 to 2 for Panel A-C) shown in white. A positive sign indicates increased 
outgassing of CO 2 to the atmosphere. 



Fig. 5. Global flux anomaly as a function of time-lag for 5 El Nino 
events from 1997-2009. El Nino events (colored lines) indexed 
by year (Panel A). Lags (/mei) from 0 to 6 months shown; flux 
anomaly always lags MEI. Flux anomalies (PgCyr -1 ) are Azotic 
(Panel A), Afl re (Panel B), A to t a i (Panel C), and Jena inversion 
(Panel D). The 1997/1998 event (blue line) cannot be calculated 
with a lag of 5 or greater as it began in April and a time delay of 
5 or more months was not observed. A positive sign indicates in- 
creased outgassing of CO 2 to the atmosphere. 


This lack of consistent response in space or time occurred 
despite a significant correlation between MEI and monthly 
terrestrial carbon sink anomalies (Fig. 8) across the full 13- 
year data record (cf. Patra et al., 2005b). This was the case 
for both the Jena inversion (r — 0.41, p < 0.001) and A[ n - C 
(r — 0.30, p < 0.001). No significant correlation was found 
using either biotic or total flux anomaly from the bottom-up 
method. Allowing for a lagged response (/mei = lug with 
maximal correlation) altered these values for the Jena inver- 
sion (/ M ei = 3, r = 0.51, p < 0.05), Afi re (/mei = 4, r = 0.47, 
p < 0.05), and Azotic (/mei = 6, r — 0.25, p < 0.05). No 
change, and no significant lag, was found for A to tai- 

In the spatial domain MEI and water deficit exhibited 
a clear teleconnection. Large-scale teleconnective clusters 
with similarly signed correlation occurred globally but were 
sparser in Russia (Fig. 9a). Similarly, Jena inversion results 
and A to tai (Fig- 9b and c) showed significant correlations 
with MEI. While both flux responses exhibited similarities 
in the aggregate, the larger grid cell (3.75° x 5°) and in par- 
ticular the a priori correlations of the Jena results masked 
smaller-scaled features in several regions, e.g., the dipole pat- 
tern in sub-Saharan Africa. Comparing lagged to unlagged 
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Fig. 6. Patterns of Abiotic (gCm 2 yr ' ) for 5 El Nino events from 1997-2009. Inset years denote El Nino event. Values are event totals 
scaled to 1 year. Non-vegetated grid cells and I A^ioticl 5 50g C m~ 2 yr -1 shown in white. For the 2009 event 61 grid cells in temperate 
Eurasia and Indochina have values from 500 to 1 166 g C m~ 2 yr - 1 . Positive anomalies indicate increased CO 2 outgassing (decrease in sink 
strength) and vice versa. 



Fig. 7. Patterns of Afl re (gCm 2 yr 1 ) for 5 El Nino events from 1997-2009. Inset years denote El Nino event; panel insets show insular 
Southeast Asia. Values are event totals scaled to 1 year. Non-vegetated grid cells and |Afj re | < 25gCm~ 2 yr -1 shown in white. Positive 
anomalies indicate increased fire emissions and vice versa. 


teleconnection maps (not shown), the highest correlation, in 
absolute value, increased from 0.59 (cf. Fig. 9b) to 0.61 and 
the median value increased from 0.39 to 0.41. 

Overall, the temporal correlation and spatial teleconnec- 
tion were not sufficient to regularize terrestrial carbon sink 
response to El Nino. While correlations were significant for 
both the top-down and bottom-up methods and approached 
0.6 in magnitude (Fig. 9b), this only equated to ~36 % vari- 
ation explained in the best case. The amount of variation ex- 


plained by El Nino was also constrained by the 13-year anal- 
ysis period. This length effectively precluded elimination of 
longer-term and possibly confounding trends or cyclical phe- 
nomena such as changes in soil moisture limitation (Jung et 
al., 2010), global dimming and brightening (Wild, 2009), and 
El Nino type based on spatial signatures of equatorward sea 
surface temperature anomalies (Collins et al., 2010; Yeh et 
al., 2009). Thus, ENSO as a mode of global climate vari- 
ability (with or without a time-lag response) cannot constrain 
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Fig. 8. Temporal profiles of global flux anomaly and MEI from 
1997-2009. Monthly time series show MEI (red), A tota i (blue), 
and Jena (green). For clarity the time series for A tota i and Jena 
were loess smoothed with a 12-month span. Lighter symbols show 
unsmoothed anomalies. Year label is centered on June-to- August. 
Grey background indicates El Nino event. Jena inversion results 
only available from 1997-2008. 


terrestrial sink behavior as the majority of variation remained 
unaccounted for. As a result the terrestrial integrated signal 
of the 5 observed El Nino events was not consistent. 

These observed dynamics were present even within the 
context of uncertainty. Spatial patterns in Sbiotic showed more 
similarity to uncertainties in NEE sensitivity as opposed to 
reanalyzed water deficit (Fig. 10). Monthly values of erotic 
were highly left-skewed with a mean of 6.5 g C m -2 month -1 
and a maximum value of 43.1 g Cm -2 month -1 , co-located 
in a cluster of relatively high values in Manitoba, Canada. 
Higher values of ebiotic generally occurred in areas domi- 
nated by croplands including the Great Plains region of North 
America, central Russia, and China (Fig. 10c). 

Uncertainty values decreased when aggregated to longer 
time scales and larger spatial domains. Average yearly 
Sbiotic by pixel, based on summation in quadrature, was 
22.5 gC m -2 yr -1 . Furthermore, integrated globally and ag- 
gregated to annualized El Nino events, uncertainty was gen- 
erally an order of magnitude less than flux anomaly (Table 1) 
and ranged from 4% (2006/2007 event) to 40% (1997/1998 
event). In all cases, 95 % confidence intervals did not include 
0, i.e., all Abiotic. Afj re , and A to tab as well as Jena (Table 1), 
single event totals were statistically non-zero. 

Incoherent responses to El Nino were also found in pre- 
vious top-down studies. For example. Bosquet et al. (2000) 
showed both an anomalous source and sink during El Nino 
years. Similarly, Patra et al. (2005a) found only 1 of 1 1 land 
regions exhibited a correlation between ENSO activity and 
CO 2 flux anomaly that translated into at least 50% varia- 
tion explained: Tropical Africa using a 3-month lag showed 
r 2 = 0.5041 or 50.41 % variation explained. We contend that 
the ongoing debate on the dynamics between drought coin- 
cident with El Nino and CO 2 flux in the Amazon needs to be 
broadened to the global scale. 



- 0.5 0 0.5 


Fig. 9. Teleconnections with MEI. Correlations between MEI and 
water deficit anomalies (Panel A), MEI and Jena inversion results 
(Panel B), MEI and Abiotic (Panel C). Correlations were computed 
on a monthly time step for each pixel relative to MEI across the full 
analysis period (1997-2009). Non-vegetated grid cells and pixels 
without significant teleconnections (p > 0.05) shown in white. 


Comparing the tower upscaling and atmospheric inversion 
approaches we found large across-event variability in the 
spatial pattern and magnitude of carbon flux responses to El 
Nino. However, there was no overlap between the 2 estimates 
of total annualized anomaly; Rodenbeck et al. (2003) anoma- 
lies were outside 95 % confidence bounds of the bottom-up 
Atotai for all events (Table 1). Spatially, large-scale regions of 
agreement (Figs. 3 and 4) were found only in southern Africa 
and northern Europe. This general lack of correspondence 
between bottom-up and top-down approaches derives from 
underlying methodological limitations and assumptions: 

1. The top-down method quantifies total surface- 
atmosphere CO 2 exchange, while the upscaling method 
was constrained to solely address El Nino-induced 
drought effects using water deficit and changes in fire 
emissions. Thus, the top-down inversion includes other 
factors, e.g., possible light and temperature drivers 
of El Nino-induced anomalies in CO 2 exchange, as 
well as any longer-term trends associated with global 
environmental change such as a general increase in 
temperatures (Meehl et al., 2007; Trenberth et al., 2007) 
and CO 2 fertilization (Friedlingstein et al., 2006). 
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Fig. 10. Component and total uncertainty in biotic flux anomaly. 
Uncertainties are (lcr) for each grid cell across all months (1997- 
2009). Panel A: Uncertainty (g C m -2 month - 1 ) in NEE sensitivity 
to water deficit (Schwalm et al., 2010a). Panel B: Uncertainty (di- 
mensionless) in water deficit based on a 10-member ensemble of 
reanalysis products (Bosilovich et al., 2009). Panel C: Total uncer- 
tainty (Sbiotio g C m -2 month 1 ) in Abiotic from combining compo- 
nent uncertainties (Panels A and B) using Monte Carlo techniques. 
Uncertainty was time-averaged with weights from GIMMS /par 
normals (Tucker et al., 2005) to emphasize periods of greater bi- 
ological activity. Non-vegetated grid cells shown in white. Note 
different color scales. 

2. The CO 2 concentrations and flux tower networks both 
drastically undersample the tropics (cf. Scholes et al., 
2009). This is problematic as the impact of El Nino 
events on the terrestrial biosphere is assumed larger in 
the tropics than elsewhere (McPhaden et al., 2006). The 
undersampling of the tropics (Fig. 1) makes it diffi- 
cult for inversions to resolve fluxes in tropical Africa 
vs. tropical South America (Gurney et al., 2003). For 
FFUXNET a denser network was used (Fig. 1): 938 
site-years across 238 sites compared to 53 sites for the 
Jena results. The highest degree of FFUXNET rep- 
resentativeness, i.e., the largest amount of site-years 
used in NEE sensitivity derivation, was in the mid- 
latitudes of the Northern Hemisphere (Fig. 11). In the 
tropics, where NEE sensitivities in evergreen broadleaf 
forests reached ~12gCm -2 month -1 (Schwalm et al., 
2010a) and water deficit anomalies ranged from -2.85a 
to +3. 28a, only ~5% of all site-years were available: 
35 site-years in evergreen broadleaf forests and 12 site- 
years in non-forested types. 


!/) 



Fig. 11. Representativeness of FLUXNET. Values are from a 
mapped cross-tabulation of site-years by Koppen-Geiger climate 
(Peel et al., 2007) and IGBP land cover class (Loveland et al., 
2001). The FLUXNET compilation used in NEE sensitivity estima- 
tion contained 938 site-years across 238 sites collected from 1991 
to 2006. Non-vegetated and non-sampled areas shown in white. 

3. Different input datasets introduce uncertainty in both 
methods. MERRA is a reanalyzed product. While 
~4 million observations are assimilated each 6-hr pe- 
riod (http://gmao.gsfc.nasa.gov/merra/) the gridded re- 
analyzed fields were ultimately generated using a land 
surface model (LSM). MERRA used the Catchment 
LSM (Koster et al., 2000), embedded in the Goddard 
Earth Observing System Model, Version 5 (GEOS-5). 
Although Catchment LSM has been shown to perform 
well in model evaluation examinations (e.g., Dirmeyer 
et al., 2006) model errors could still propagate through 
the bottom-up approach. For the Jena inversion results 
uncertainty was in the transport model and meteorologi- 
cal forcing data, in addition to the uncertainty from data 
sparseness and a priori assumptions. Thus, we quan- 
tified the uncertainty associated with reanalyzed water 
deficit and propagated this through to flux anomalies 
and showed that our findings were robust after discount- 
ing this uncertainty as well as the uncertainty associated 
with the Jena results (Table 1). 

4. Sub-grid variability for the Jena product was larger than 
for the bottom-up approach based on their respective 
resolutions: 3.75° x 5o for Jena and 1° x 1° for the 
bottom-up approach. More important than the pixel 
size is the larger a priori correlation length which forced 
the inversion to be smooth on at least a 1200 km scale 
(Rodenbeck et al., 2003). This, coupled with known 
issues of spatial resolution in inversion-constrained es- 
timates (e.g., Gurney et al., 2003), biased the spatial 
intercomparisons. Consistency did not improve, both 
across events and methods, even when both methods 
were coarsened beyond the correlation length (Fig. 4). 

Our results clearly show that El Nino events do not produce 
a consistent CO 2 response from the terrestrial biosphere, 
even though top-down and bottom-up estimates did not show 
statistical agreement in response. High variability in total 
biospheric anomaly, both across-method and within-method. 
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underscored the general lack of consistent pattern of terres- 
trial sink response to El Nino. In addition, both bottom- 
up and top-down methods revealed that El Nino does not 
impose a coherent across-event response regarding the inte- 
grated global terrestrial response as well as its spatial pattern- 
ing. Finally, even though monthly time series of global net 
terrestrial-atmosphere CO 2 flux anomalies were significantly 
correlated with MEI for both approaches (flux upscaled and 
atmospheric inversion), this did not result in the canonical ex- 
pectation of increased CO 2 outgassing during El Nino events. 

4 Conclusions 

There is no consistent response of terrestrial carbon cycling 
to El Nino events. The significant correlation between MEI 
and global CCEflux anomalies, estimated using either ap- 
proach, did not translate into spatiotemporal consistency in 
terrestrial sink response to El Nino. Although we cannot ex- 
clude confounding longer-term trends or periodicities based 
on the time domain of this analysis, our results indicate that 
drought is insufficient to constrain the effects of El Nino 
on carbon cycling. While ENSO is the dominant mode 
of year-to-year climatic variability, this was not reflected 
in global CO 2 flux anomalies. Importantly, the bottom-up, 
data-oriented approach presented here indicates that El Nino 
events may act to either enhance or diminish terrestrial car- 
bon uptake. This is in contrast to the prevailing view that El 
Nino induces CO 2 releases by stimulating drought. Further- 
more all 5 El Nino events across the 13-year record displayed 
a unique Abiotic spatial signature with only limited coherence 
in parts of Amazonia, Australia and southern Africa. Across 
all events, fire emissions anomalies were smaller in magni- 
tude than biotic anomalies of CO 2 exchange except during 
the 1997/1998 El Nino event, i.e., the global terrestrial re- 
sponse to El Nino is driven more by the biotic component 
than by fire. Lastly, the bottom-up scaling approach detailed 
here is a biologically interpretable and efficient method to 
extend the kernel of observed FLUXNET data to large-scale 
spatially explicit patterns of carbon cycle dynamics. 
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